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Lattice QCD calculations of hadron structure have been carried out for two decades, see [1] 
for a recent review. As the up- and down-quark masses are gradually lowered toward the physical 
point (isospin symmetry is assumed), new questions have arisen on the way. The calculations have 
reached a certain degree of maturity and consensus among different collaborations down to about 
m % = 300MeV. With the goal in mind of confronting the lattice results with hadron phenomenology, 
it is useful to first discuss the quantities that are well determined experimentally. Such quantities 
include the axial charge of the nucleon, the isovector momentum fraction and the electromagnetic 
(e.m.) form factors (FFs). One may then focus on quantities where the lattice can be a discovery 
tool, such as the generalized form factors. 

I therefore start by discussing the calculation of e.m. FFs, both because they are a special case 
of the GFFs and because they represent a benchmark calculation of experimentally well-determined 
functions. The matrix elements of the e.m. current between two nucleon states are parametrized by 
the Dirac and Pauli form factors (q = P' — P, Q 2 = —q 2 ), 

(P',A'\^\P,A)=U(P',A')r^q 2 )U(P,A), (q 2 ) = f F X (Q 2 ) + io^ ' ^F 2 {Q 2 ) . (1) 

The result of a calculation at m % = 300MeV by the LHP collaboration is displayed in Fig. (1, top 
row). It uses 2+1 flavors of domain- wall fermions at a lattice spacing of a = 0.084fm. Several 
other collaborations have obtained the FFs at similar pion masses [2]. A common point among 
these calculations is that the dipole form fits the data well, but the radii extracted from the fits are 
significantly smaller than the phenomenological radii. Secondly, chiral effective theory predicts 
a strong pion mass dependence of the Dirac and Pauli radii, but the lattice data exhibits a much 
milder dependence. It may be that the ra^-range of applicability of the chiral formulae is much 
smaller than 300MeV. Finite-volume effects need to be investigated in more detail too. 

The approximate dipole behavior of the phenomenological FF F\{Q 2 ) = l+ Q2^ llGeW 2 can 
be understood as being due to the contribution of two nearby vector meson poles with opposite 
residua [3]. Figure (1, bottom row) illustrates the fact that in Nature, the Dirac dipole mass lies 
within 10% of the p meson mass. In QCD at larger m u ^ however, lattice results show that the Dirac 
dipole mass is significantly larger than the vector meson mass in the same theory. We note that, 
since both the Dirac and Pauli radii diverge in the chiral limit, the latter must become very large 
compared to the p Compton wavelength approaching the chiral limit. 

The numerical evaluation of Wick-connected and -disconnected diagrams that contribute to 
three -point functions (Fig. (1), bottom row) proceeds in a rather different way. Most collaborations 
have focused on the more tractable connected contributions, thus restricting themselves to isovector 
quantities. Consider in this respect the flavor structure of the e.m. FFs of the proton 1 , 

77 _ 2 pu _ 1 pd _ 1 i7.v _ 1 pu-d , 1 pu+d-ls _ 1 pu-d i 1 pu+d , i i pu+d _ j ps \ 
1 1,2 — 3^ 3^ 3^ — i ^6 — 2 -"conn ' 6 conn ' 6 v disc ^discJ- 

Figure (1, middle row), displaying the strangeness FFs [4], suggests that disconnected diagrams 
contribute less than 0.01 to F\, which is negligible at low Q 2 compared to the other uncertainties. 

I now briefly describe the connection between the generalized parton distributions (GPDs), 
which govern Deeply Virtual Compton Scattering, and the GFFs that can be computed on the 

'NB. The normalization is such that F"(0) = 2, and the Sachs FFs are given by Ge = F\ — Gm = F\ +F2. 
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Figure 1: Top row: the e.m. FFs in an Nf = 2 + 1 calculation using domain-wall fermions at a lattice spacing 
0.084fm (0.1 14fm in the 'coarse' case). The right plot displays the Pauli FF calculated at m K = 300MeV [5]. 
Middle row : calculation of the strangeness FF of the nucleon by the %QCD collaboration [4] (Nf = 2 + 1, 
a = 0.12fm, m K > 600MeV). Bottom row : connected and disconnected Wick contractions contributing to 
hadronic matrix elements (fig. from [4]). Right, the nucleon Dirac dipole mass [6]. 



lattice. The GPDs parametrize the matrix elements of the operators 



ig J_'^j 2 dan-A(an 



*( T )- 



(2) 



Here n is a light cone vector, fi = , and in the following x is the momentum fraction, A = P' — P, 
t = A 2 and t, = —n ■ A/2. Different F matrices lead to different GPDs, 

T = fi^(H,E), r = fa<^(H,E), r = n^ j ^(H T ,E T ) and (H T ,E T ). (3) 

For r = fi, the 'unpolarized' case, the parametrization at a renormalization scale \i takes the form 

(P\A'\^(x^)\P,A) = ((m(^,t,^ + ^((ian)E(x^,t^), (4) 
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where ((r)) = U (/ y , A')TU (P,A) and U (P,A) is a spinor that solves the free Dirac equation. GPDs 
determine the distribution of partons and their helicities in impact parameter space [7], 



q(x,b ± ) = I^e- iA ^H(x^=0-A 2 ± ), Aq(x,b ± ) = J^e- iA ^H(x^ = 0,-A 2 ). 



Mellin-moments of the GPDs, H n (£,t) = j\Axx n - l H(x,t,,t), are given by polynomials in 
the longitudinal momentum transfer E, , 



The coefficients of these polynomials are the GFFs calculated on the lattice. In the n = 1 case, 
A\q and Bio coincide with F\ and F2. In the n = 2 unpolarized case, the operator is ff^ lfl2 (x) = 
yiY^D^ + Y^D^q and, with P = {P + P')/2, 



(^|^ lM2 (^)|P)=P {w ((7 M2} ))A 2 o(0 + ^ 1 ((a M2}a ))A a B 20 (0 + ^A {Ml A / , 2} C 2 o(0. (6) 



The forward-limit of A 2 o and B20 determine the momentum and spin decomposition of the nu- 
cleon in terms of quark and gluon contributions [8]. In practice, only the first few moments n < 3 
are accessible on the lattice. The impact-parameter dependence bj^ can be obtained more easily, 
currently (0.4GeV) 2 < A 2 < (2.0GeV) 2 . The lower end of this interval is limited by momentum 
quantization in a finite box, P = 2 j^, k € Z. It can potentially be overcome with twisted boundary 
conditions, whereby P = £ + ^ [9]. The upper end is limited by the inverse lattice spacing. 

Figure (2, left panel) displays the n = 1,2 and 3 GFFs of the nucleon (isovector, unpolarized 
case). This particular calculation is performed in a (3.5fm) 3 box at a pion mass of 355MeV. For the 
purpose of the comparison, the FFs have been normalized by their forward value. It is clearly seen 
that the FFs become flatter as n increases. This translates into decreasing transverse radii, (right 
panel). The latter have been obtained by performing dipole fits to the FF data with \t\ < 0.5GeV 2 . 
The simulations are performed in boxes of size (2.5fm) 3 , with in addition a (3.5fm) 3 simulation 
at m n = 355MeV. In the n = 1 case, a gradual increase in the (Dirac) radius is observed as the 
pion mass is lowered toward its physical value. For n>2, this growth in the transverse size of the 
nucleon is not observed in the range of explored pion masses. This is qualitatively consistent with 
the expectation from chiral effective theory, which predicts a logarithmic divergence in the n = 1 
case, but finite values for the radii corresponding to higher n. However, at m % = 355MeV one also 
observes a statistically significant difference between the radii extracted from volumes (2.5fm) 3 
and (3.5fm) 3 . Since the finite-size effect is expected to increase when the pion mass decreases, it 
appears necessary to check the observed trends by using large-volume simulations. 

In conclusion, GPDs play a central role in describing quantitatively the three-dimensional 
structure of the nucleon. Lattice QCD is ideally suited to determine the GFFs associated with their 
low Mellin-moments. A handful of lattice collaborations have published results for the nucleon 
structure for pion masses down to about 300MeV. A general observation is that many observables 
show a very mild pion mass dependence. The non-analytic pion-mass dependence predicted by 
chiral effective theory is not yet visible in the data; presumably it sets in at a lower pion mass. The 
next step will therefore be to go down to 200MeV [10] with controlled uncertainties, in particular 
a solid understanding of the finite-volume effects in the simulations is crucial. 



H"=\^,t) =A 10 (0, 

H"= 2 (^,t) =A 20 (0 + (2£) 2 C 2 o(0 

H"='^,t) =A 30 (t) + (2^) 2 A 32 (t) 



E n =\^t) =B l0 (t), 

E"= 2 ^,t) =S 20 (0-(2^) 2 C 20 (0, 

E"= 3 ^,t) =B 30 (t) + (2t;) 2 B 32 (t), ... 



(5) 
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Figure 2: Left: Unpolarized isovector GFFs. Right: the corresponding radii [6]. 

The formalism to study transverse-momentum dependent PDFs on the lattice, relevant to semi- 
inclusive DIS, has recently been studied [11]. One important aspect of the calculation is the renor- 
malization of non-local operators of the type q(x)rU x ,oq(0), where U x fl is a Wilson line. A second 
essential issue is the connection between the matrix elements of this operator with U x q along the 
light-cone and with U in the Euclidean domain. I expect to see further developments in this field 
in the near future. 
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